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ABSTRACT 

An analysis is presented of the power spectrum of X-ray variability of the bright Seyfert 1 
galaxy Mrk 766 as observed by XMM-Newton. Over the 0.2-10 keV energy range the power 
spectral density (PSD) is well-represented by a power-law with a slope of aiow ~ 1 at low 
frequencies, breaking to a slope of ckhi = 2.8^04 at a frequency /bi ~ 5 x 10"'* Hz. As has 
been noted before this broken power-law PSD shape is similar to that observed in the Galactic 
black hole candidate Cygnus X-1. If it is assumed that Mrk 766 shows a power spectrum 
similar in form to that of Cyg X-1, and that the break timescale scales linearly with black 
hole mass, then the mass of the black hole in Mrk 766 is inferred to be < 5 x 10^ M0. 
This rather low mass would mean Mrk 766 radiates above the Eddington limit. The coherence 
between different energy bands is significantly below unity implying that variations in the 
different energy bands are rather poorly correlated. The low coherence can be explained in the 
framework of standard Comptonisation models if the properties of the Comptonising medium 
are rapidly variable or if there are several distinct emission sites. 

Key words: galaxies: active - galaxies: Seyfert: general - galaxies: individual: Mrk 766 - 
X-ray: galaxies 



1 INTRODUCTION 

X-ray variability appears to be ubiquitous in Active Galactic Nuclei 
(AGN). The rapid and seemingly random fluctuations in the X-ray 
luminosity of Seyfert galaxies provided early support for the stan- 
dard black hole/accretion disc model (e.g. Rees 1984) by implying 
compact emission regions and high luminosity densities (Fabian 
1979; M'^Hardy 1989). 

The long (~ 3 day), uninterrupted observations possible with 
EXOSAT allowed the power spectral density (PSD) of Seyfert 1 
galaxies to be measured for the first time (Lawrence & Papadakis 
1993; Green, M'^Hardy & Lehto 1993). The PSD represents the 
average of the (squared) amplitude of variations as a function of 
temporal frequency (see van der Klis 1989). The £X05Ar observa- 
tions showed that the PSDs of Seyfert Is above ~ 10^^ Hz could 
be approximated by a power-law {V{f) oc f^" where is 
the power at frequency / and a is the PSD slope) rising steeply at 
low frequencies with a slope a > 1 (so-called "red noise" spec- 
tra; Press 1978). It was noted early on (e.g. Lawrence et al. 1987; 
M'^Hardy 1989) that this red noise variability of Seyferts is similar 
to that observed in Galactic Black Hole Candidates (GBHCs; van 
der Klis 1995), perhaps suggesting that the same physical processes 
operate in these sources that differ in black hole mass by factors of 
> 10^ 

The steep slopes of the £X05ArPSDs mean that at lower fre- 
quencies the PSD must flatten in order for the integrated power to 
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converge. Indeed, breaks in the PSDs of Seyfert galaxies have re- 
cently been detected (Edelson & Nandra 1999; Uttley, M'^Hardy 
& Papadakis 2002; Vaughan, Fabian & Nandra 2003; Markowitz 
et al. 2003). The position of these breaks represent "characteristic 
timescales" in the aperiodic variability of Seyfert Is, and interest- 
ingly they appear to scale linearly with the mass of the central black 
hole (Markowitz et al. 2003). The PSDs of GBHCs usually show 
at least two breaks (Nowak et al. 1999; Lin et al. 2000), and it is 
becoming increasingly apparent that the PSDs of Seyfert galaxies 
(and their timing properties in general) are very similar to those of 
GBHCs except shifted to longer timescales according to the black 
hole mass (Vaughan et al. 2003; Markowitz et al. 2003). 

This letter presents an analysis of the X-ray continuum vari- 
ability of the bright Seyfert 1 galaxy Mrk 766 (aka. NGC 4253; 
z — 0.012929) using a long XMM-Newton observation. XMM- 
Newton offers the possibility to improve upon the high frequency 
PSDs provided by EXOSAT, having as it does a similarly long or- 
bit but greatly increased throughput and a wider bandpass. The rest 
of this paper is organised as follows. Section|2|describes the basic 
data reduction procedures, Section|3|details the results of the PSD 
analysis and Section|4|describes the cross spectral results. Finally 
the implications of these are discussed in section|5| 

2 DATA REDUCTION 

Mrk 766 was observed by XMM-Newton (lansen et al. 2001) over 
the period 2001 May 20 - 2001 May 21 (rev. 265), during which 
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Figure 1. Full-band (0.2-10.0 keV) EPIC pn light curve in 100 s bins. 
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all instraments were operating nominally. The present analysis is 
restricted to the data from the European Photon Imaging Cameras 
(EPIC). Other aspects of this observation are discussed in Mason et 
al. (2002). The EPIC pn camera (Struder et al. 2001) was operated 
in small-window mode, as was the M0S2 camera (Turner et al. 

2001) , and MOSl was operated in timing mode. Here only data 
from pn are used as these have the highest signal-to-noise ratio and 
are free from photon pile-up. 

Extraction of science products from the Observation Data 
Files (ODFs) followed standard procedures using the XMM- 
Newton Science Analysis System v5.3.3 (SAS). The EPIC data 
were processed using the standard SAS processing chains, source 
data were extracted from a circular region of radius 35 arcsec from 
the processed pn image. Only events corresponding to patterns 0- 
4 (single and double pixel events) were used for the pn analy- 
sis. Background events were extracted from regions in the small 
window least effected by source photons. These showed the back- 
ground to be relatively low and stable throughout the first 105 ksec 
of the observation. During the final ~ 15 ksec the background rate 
increased dramatically. The much stronger (and highly variable) 
background level during this part of the observation contributed a 
significant amount of non-source variability, and so only the first 
105 ksec of the light curve are used for the variability analysis pre- 
sented here. 

Light curves were extracted from the EPIC pn data in four dif- 
ferent energy bands: 0.2-10.0 keV (full band), 0.2-0.7 keV (soft 
band), 0.7-2.0 keV (medium band) and 2.0-10.0 keV (hard band). 
These were corrected for telemetry drop outs (less than 1 per cent 
of the total time), background subtracted and binned to 100 s time 
resolution. The errors on the light curves were calculated by prop- 
agating the Poisson noise. The light curves were not corrected for 
the ~ 71 per cent "live time" of the pn camera (Struder et al. 2001), 
which is only a scaling factor. The full band light curve is shown 
in Fig. □ The average source (background subtracted) and back- 
ground count rates are shown in Table Q along with the fractional 
excess rms variability amplitude of the source (Fvar; Edelson et al. 

2002) in each energy band. 



3 POWER SPECTRAL ANALYSIS 

The analysis presented here followed very closely that applied by 
Vaughan et al. (2003) to XMM-Newton light curves of the bright 
Seyfert 1 galaxy MCG-6-30-15. The periodogram of the con- 
tiguous (100 s binned) light curve was calculated using the stan- 



Table 1. Basic properties of the light curves. Mean source and background 
count rates and fractional excess rms variability amplitudes (i^var) in each 
energy band. 



Band 


source (ct s ^ ) 


background (ct s ^ ) 


Fvar (%) 


Full 


19.9 


0.3 


17.6 ±0.4 


Soft 


11.2 


0.2 


17.8 ±0.4 


Medium 


6.68 


0.07 


19.0 ±0.4 


Hard 


2.03 


0.05 


17.8 ± 0.4 



dard Discrete Fourier Transform (DFT). The periodogram was then 
binned using the method of Papadakis & Lawrence (1993), using 
A'^ = 15 periodogram points per bin, to produce a consistent PSD 
estimate with Gaussian error bars The periodogram was normalised 
to "(rms/mean)^ Hz~^" units as defined by van der Klis (1997). 
The binned periodogram is shown in Fig. |2| and appears steep at 
low frequencies (due to the red noise variability of the source) and 
flattens at high frequencies due to the variability introduced by the 
Poisson noise (which has a flat, "white" power spectrum). 

The periodogram was fitted using the Monte Carlo procedure 
described in Vaughan et al. (2003) (similar to that discussed by Utt- 
ley et al. 2002) in order to derive parameters of the PSD of Mrk 766. 
For each trial PSD model 500 simulated light curves were gener- 
ated (using the algorithm of Timmer & Konig 1995). The simulated 
light curves were generated from a PSD model that extended down 
to much lower frequencies than the measured periodogram, and 
each simulated light curve was a factor ~ 50 longer than the ob- 
served light curve. This allows for variability power on timescales 
longer than the observed light curve and thereby accounts for any 
effect this might have on the observed periodogram ("red noise 
leak" - see Uttley et al. 2002 and Vaughan et al. 2003). A sec- 
tion of each long simulated light curve was re-sampled to match 
the sampling of the observed light curve. Periodograms were then 
calculated for each of the 500 simulated datasets (and binned ex- 
actly as the real data were) and the 500 binned periodograms aver- 
aged to produce an average periodogram. This average represents 
the PSD model after being distorted by the light curve sampling 
("folded" through the response of the observation). The constant 
Poisson noise component is then added (at the level expected for 
Poisson noise) and the folded model is compared to the data by 
measuring the of the fit. Further details of the procedure are 
given in Vaughan et al. (2003). 

Fitting the data with a power-law PSD model gave a poor fit 
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Figure 2. Top panel: binned periodogram of the full band (0.2-10 keV) 
light curve (data points) and best-fitting single power-law model (soUd line). 
The dotted line marks the expected Poisson noise level. Bottom panel: resid- 
uals from the power-law fit. Inset panel: variation in with power-law 
slope. 
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Figure 3. Top panel: binned periodogram of the full band (0.2-10 keV) 
light curve (data points) and best-fitting broken power-law model (solid 
line). Bottom panel: residuals from the broken power-law fit. Inset panel: 
68.3, 90 and 99 per cent confidence contours for the model parameters 
(break frequency and power-law slope above the break). 



to the data (x^ = 80.8 for 33 de grees for freedom dof, rejected 
at > 99.9 per cent confidence) with a best-fitting power-law slope 
of a = 1.95. The steep slope seen from these XMM-Newton data 
requires that the PSD flatten at lower frequencies (as observed in 
other Seyfert 1 galaxies; Uttley et al. 2002; Vaughan et al. 2003; 
Markowitz et al. 2003) for the total integrated power to be finite. 
The poor fit of the single power-law model further suggests that the 
power spectral break may occur within (or close to) the frequency 
range probed by this XMM-Newton observation. 

In order to test this hypothesis a broken power-law PSD model 
was fitted to the data. The PSD was assumed to break from a slope 
of aiow = 1 at low frequencies to a steeper slope. The free parame- 
ters of the model were the PSD normalisation, the break frequency 
(/br) and the high frequency slope (ahi). As shown in Fig.|3] this 
model provided a reasonable fit to the data (x^ = 43.2/32 dof 
with a rejection probability of 91.5 per cent), much better than 
the unbroken power-law (Ax^ ~ 37.6 improvement for one ad- 
ditional parameter). The slightly high rejection probability may be 
indicating that a broken power-law is only an approximation to the 
true underlying PSD. The best-fitting parameters are as follows: 
Qhi = 2.81q'4 ^id break frequency of /br ~ x 10^* Hz. 
(The errors represent 90 per cent confidence limits calculated using 
a Ax^ = 2.7 criterion.) 

The PSD slope below the break is not well-constrained with 
these XMM-Newton data. An alternative model, in which the 
power-law breaks from aiow = to a steeper slope was also 
compared to the data. The best-fitting high frequency slope was 
Qhi = 2.4l°;^ and the break frequency was 1.6 ± 0.8 x 10"'' Hz, 
giving = 48.1/32 dof (rejected at 96.6 per cent), slightly 
worse than the model assuming a break from aiow = 1- Thus it 
is not possible to constrain in detail the shape of break in Mrk 766, 
but the broken power-law models did provide a significant improve- 
ment over the single power-law model. As a low frequency slope 
of 1 fitted the data slightly better, and seems more likely in other 
Seyfert galaxies (Uttley et al. 2002; Markowitz et al. 2003), this 
model is used throughout the rest of this paper to describe the PSD 
of Mrk 766. 

Energy dependence of the high frequency PSD, with the PSD 
slope being flatter at higher energies, has previously been seen in 



Table 2. Results of fitting the periodograms derived from the four energy 
bands with the broken power-law model (breaking from a slope of Qiow = 
1 at low frequencies to Ohi)- Errors on the model parameters correspond to a 
90 per cent confidence level for one interesting parameter (i.e. Ax^ = 2.7). 



Band 


Slope 

"hi 


/br 

(10-"' Hz) 


x'/dof 


full 
soft 

medium 
hard 


^■*'-0.4 

2.6 ±0.4 

2.7 ±0.5 
2.3 ±0.4 


+ 3 

5±3 
5±3 


43.5/32 
44.9/32 
51.1/32 
44.0/32 



the Seyfert 1 galaxies NGC 7469 (Nandra & Papadakis 2001), 
MCG-6-30-15 (Vaughan et al. 2003) and NGC 4051 (Papadakis 
& Lawrence 1995; M'^Hardy et al. 2003). In order to test for any 
such variation in the PSD as a function of energy the above analysis 
was repeated on the light curves extracted in the three energy sub- 
bands. The results are shown in Table|2| The high frequency slopes 
are formally consistent with one another although the hardest band 
does show the flattest best-fitting slope (this is also the case when 
the data were fitted with the PSD model assuming aiow = below 
the break). This hints that Mrk 766 may show the energy depen- 
dence seen in other Seyferts, and in particular adds to the evidence 
provided by XMM-Newton (see Vaughan et al. 2003 and M'^Hardy 
et al. 2003) that the PSDs of Seyferts are energy-dependent even 
at soft X-ray energies (the energy dependence in NGC 7469 could 
only be examined above 2 keV, the same is true for most GBHCs). 



4 CROSS SPECTRAL ANALYSIS 

The cross spectrum is related to the PSD (sometimes known as the 
auto spectrum). It offers a comparison of two simultaneous light 
curves (e.g. from different energy bands) as a function of temporal 
frequency. The cross spectrum can be used to derive the coherence 
(7^(/)) of two light curves (a measure of how well correlated the 
variations in the two light curves are) and time delays between the 
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Figure 4. Coherence function between light curves in the thi'ee energy 
bands. The coherence (at a given frequency) decreases as the separation 
between energy bands increases. At frequencies above 10" Hz the Poisson 
noise power becomes comparable to the source power and the coherence 
becomes difficult to measure (Vaughan & Nowak 1997). 



light curves (see Vaughan & Nowak 1997 and Nowak at al. 1999 for 
an explanation of these properties in terms of GBHC observations 
and Papadakis, Nandra & Kazanas 2001 and Vaughan et al. 2003 
for an application to Seyfert galaxies). 

The coherence of Mrk 766 was measured using the recipe 
given by Vaughan & Nowak (1997). The periodograms of the two 
light curves and their cross periodogram were computed and binned 
by TV = 15 consecutive frequencies per bin The coherence and its 
uncertainty were then calculated using equation 8 of Vaughan & 
Nowak (1997). (Vaughan et al. 2003 have shown that uncertainty 
estimates calculated in this way are accurate for data very similar 
to those examined here.) 

The resulting coherence functions (Fig. |4j show that, partic- 
ularly between the soft and hard bands, the coherence is signifi- 
cantly below unity at all frequencies. This means that either the 
relation between soft and hard bands is non-linear or that a fraction 
(1 — 7'^(/)) of the hard band variance cannot be accounted for by 
variations in the soft band. The apparently poor (linear) correlation 
between variations in different bands is confirmed by an analysis 
of the cross correlation function, estimated using the Discrete Cor- 
relation Function (DCF; Edelson & Krolik 1988). The DCF for the 
soft versus hard band light curves reaches a peak (at zero time de- 
lay) of only 0.65. It should be noted that the reduced coherence is 
not an artifact of Poisson noise in the data as the recipe of Vaughan 
& Nowak (1997) correctly accounts for this. The low coherence 
means it is not possible to reliably search for time delays between 
the different bands. 

The hardness ratios were calculated to further examine the 
relationship between variations in different energy sub-bands 
(Fig. HJ. The hardness ratios are significantly variable but show 
little correlation with the total source flux. There is a slight trend 
for the 0.7-2.0/2.0-10.0 keV ratio to soften as the source brightens 
(as has been seen in other Seyfert Is) but the ratios show a great 
deal of scatter (much larger than expected from the size of the error 
bars) not obviously associated with changes in the full band flux. 
This complements the coherence analysis and shows the the rela- 
tionships between energy sub-bands are complex. 
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Figure 5. Hardness ratios (in 400 s time bins). Top panel: the full band light 
curve Lower panels: (left) hardness ratios plotted against time and (right) 
against full band count rate. 



5 DISCUSSION 

The shape of the high frequency power spectrum of Mrk 766 is con- 
sistent with a power-law that breaks from a flat slope (aiow ~ 1) 
to a steep slope (ahi ^ 2) at a frequency ~ 5 x 10~* Hz. There 
is no evidence for discrete features such as quasi-periodic oscilla- 
tions (previously claimed by Boiler et al. 2001, but see also Ben- 
lloch et al. 2001). Similar broken power-law PSDs have been ob- 
served in other Seyfert Is (Uttley et al. 2002; Vaughan et al. 2003; 
Markowitz et al. 2003; M'^Hardy et al. 2003). This shape for the 
high frequency PSD is similar to that of GBHCs (Nowak et al. 
1999; Lin et al. 2000). Not only is the shape the same but the am- 
plitudes are similar. The power (in f'P{f) units) measured at the 
break frequency is ~ 0.01 in the 2-10 keV band, which compares 
with the power at the high frequency break in Cyg X-1 of ^ 0.02 
(Belloni & Hasinger 1990). This strongly supports the idea that the 
variability mechanisms, and perhaps the X-ray emission processes 
in general, are the same in AGN and GBHCs. 

In the best studied GBHC, Cygnus X-1 in its low/hard state, 
the PSD breaks from a ~ 1 to a steeper slope at ~ 3 Hz (Bel- 
loni & Hasinger 1990; Nowak et al. 1999). In its high/soft state 
the position of the high frequency break increases to ~ 10 Hz 
(Churazov & Gilfanov & Revnivtsev 2001). It is not clear which 
of these states is closer to the behaviour of Seyferts. The corre- 
sponding break frequency in Mrk 766 is a factor of ~ 6 x 10'' or 
~ 2 X 10* times higher than that seen in Cyg X-1, depending on 
whether the low/hard or high/soft state is used for comparison. As- 
suming a 10 M0mass black hole in Cyg X-1 (Hererro et al. 1995), 
and further assuming that the characteristic break timescale scales 
linearly with black hole mass (see Markowitz et al. 2003) yields 
estimates for the mass of the black hole in Mrk 766 of ~ 6 x 10* 
or ~ 2 X 10^ M0 (comparing with low/hard and high/soft state 
PSDs respectively). These estimates are unusually low compared 
to the masses expected for Seyfert nuclei (e.g. Wandel, Peterson 

6 Malkan 1999). Indeed, Wandel (2002) estimated the mass of 
Mrk 766 to be ~ 10^ M0 based on its luminosity and optical 
broad-line properties. 

The highest mass estimate consistent with a linear scaling of 
timescales, obtained using the 90 per cent lower limit on the break 
frequency in Mrk 766 and scaling by the break seen in the high/soft 
state in Cyg X-1, is ~ 5 x 10^ M©. The (unabsorbed 0.2-10 keV) 
X-ray luminosity of Mrk 766 is ~ 3 x 10*^ erg s~^ (assuming 
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Ho = 70 km Mpc~^ and go = 0.5). Taking a fairly conserva- 
tive estimate of the bolometric luminosity of 10** erg s~^ and the 
highest black hole mass estimate suggests Mrk 766 is radiating at 
~ 1.5I/Edd. Of course, this depends on there being a linear scal- 
ing between the break timescales in Cyg X-1 and those of Seyfert 
Is. If this is the case the implication seems to be that Mrk 766 is 
accreting above the Eddington limit. 

A further point is that if the slope of the low frequency PSD is 
significantly different from the assumed aiow = 1 then the break 
frequency could be lower (fitting a model assuming qiow = gave 
/br = 1.6 ± 0.8 X 10"* Hz). Comparing this to the low-frequency 
break in Cyg X- 1 (the break from a ~ to a ~ 1 in the low/hard 
state) gives an even lower black hole mass estimate (see the discus- 
sion in Uttley et al. 2002). However, the fact that the high frequency 
PSD is steep (« 2.3 in the hardest band) suggests a slope of ~ 1 
below the break is more likely, by comparison with the PSD of Cyg 
X-1. 

It is perhaps worth noting that Mrk 766 has previously been 
classified a Narrow-line Seyfert 1 (NLSl), a class of AGN of- 
ten suggested to be accreting close to the Eddington limit (Boiler, 
Brandt & Fink 1996). The PSDs of other NLSls also seem to show 
break frequencies somewhat higher than found in "normal" Seyfert 
Is (Ark 564: Pounds et al. 2001; Papadakis et al. 2002; NGC 4051: 
M'^Hardy et al. 2003). This is consistent with the idea that NLSls 
harbour lower mass black holes accreting at a higher rate (relative 
to Eddington) compared to the normal Seyfert population. 

It remains plausible that the mass-timescale relation is not per- 
fectly linear. The most important physical timescales associated 
with accretion flows (e.g. section 5.8 of Frank, King & Raine 1985) 
derive from a distance and a speed over which some physical pro- 
cess operates (such as viscosity). If the salient distance (such as the 
radius of the innermost edge of the accretion disc) is approximately 
the same in units of the gravitational radius (rg = GM/c?) be- 
tween Cyg X-1 and AGN then the timescales should (to first order) 
scale linearly with black hole mass. However, if the relevant dis- 
tance is not constant then this linear scaling can break down. For 
example, the innermost stable circular orbit around a black hole de- 
pends on its dimensionless spin parameter (a/M), with a spinning 
black hole allowing a steady accretion disc to extend to smaller 
radii (in gravitational units). If Mrk 766 has a higher spin param- 
eter than Cyg X-1 one would perhaps expect the accretion disc to 
extend to smaller rg and the physical timescales to be somewhat 
smaller than the linear expectation. Such a situation would allow 
the black hole mass of Mrk 766 to be higher by a factor of a few 
(compared to the linear assumption) and thereby allow for a sub- 
Eddington accretion rate. An independent estimate of the mass of 
Mrk 766 (from e.g. reverberation mapping or gas kinematic stud- 
ies) would help clarify this issue. 

The variations in different energy bands show a coherence 
well below unity at all frequencies examined. A loss of coher- 
ence at high frequencies has also been observed in MCG-6-30-15 
(Vaughan et al. 2003), but at the lowest frequencies probed the co- 
herence was much higher than observed in Mrk 766. On longer 
timescales the variations in NGC 7469 were consistent with unity 
coherence (Papadakis et al. 2001). In Cyg X-1 the coherence is 
known to be high (near unity) over a wide range in frequencies 
but falls significantly below unity at high frequencies (Nowak et 
al. 1999). The reduced coherence means the relationship between 
variations in different energy bands cannot be described by a single 
(linear) transfer function. 

If the X-rays are produced by inverse-Compton scattering in 
a hot corona, as is often thought (e.g. Haardt & Maraschi 1991), 



then the low coherence implies rapid changes in the coronal prop- 
erties and/or multiple independent emission sites. As discussed by 
Nowak & Vaughan (1996), a static Comptonising corona should 
produce unity coherence. Rapid changes in the properties (e.g. 
size, temperature, optical depth) of the X-ray emitting region could 
cause changes in the transfer function (relating soft to hard flux) on 
the timescales probed and so reduce the coherence. Alternatively 
the low coherence may mean there is more than one transfer func- 
tion operating simultaneously, which could be physically realised 
if the X-ray emission is produced at several independent sites. 

XMM-Newton has opened the door on high frequency timing 
studies of AGN, a door that has remained largely closed since the 
end of the EXOSAT mission. Future work must establish whether 
the high frequency power and cross spectral properties observed 
in the handful of Seyfert Is with XMM-Newton long-looks (see 
also Vaughan et al. 2003 and M'^Hardy et al. 2003) are universal 
or whether different objects show different characteristics (perhaps 
as GBHCs show distinct states). 
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